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HIGHLIGHTS 


►  A  1  kWe  thermoelectric  stack  is  considered  for  geothermal  heat  conversion. 

►  It  comprises  an  array  of  Bi-Te  modules  between  counter-flow  water  heat  exchangers. 

►  A  coupled  conduction-convection  heat  transfer  model  was  formulated. 

►  Optimization  is  accomplished  for  either  maximum  efficiency  (rj)  or  minimum  volume  (V). 

►  For  hot  water  inlet/outlet  temperatures  (413  K/393  K),  tj  =  4.2%  and  V=  0.0021  m3. 


ARTICLE 


INFO 


ABSTRACT 


Article  history: 

Received  27  January  2012 
Received  in  revised  form  8  May  2012 
Accepted  10  May  2012 
Available  online  28  June  2012 


Keywords: 

Thermoelectric 

Stack 

Geothermal 


A  thermoelectric  stack  comprising  an  array  of  Bi-Te  based  thermoelectric  converter  (TEC)  modules 
is  considered  for  geothermal  heat  conversion.  Each  TEC  module  consists  of  127  (Bi0.2Sb0.8)2 
Te3/Bi2(Teo.96Seo.o4)3  p/n- type  thermoelement  pairs,  fastened  by  30  x  30  mm2  A1203  plates.  The  thermo¬ 
element  pairs  have  leg  cross-section  of  1.05  x  1.05  mm2,  a  figure-of-merit  equal  to  1,  and  a  theoretical 
heat-to-electricity  conversion  efficiency  of  ~5%  when  the  module  is  operated  at  a  temperature  difference 
of  200  K.  A  temperature  gradient  across  the  thermoelement  legs  within  an  array  is  imposed  via  a  Cu  par¬ 
allel-plate  heat  exchanger  adhering  to  the  A1203  plates  and  operating  hot  and  cold  water  in  counter-flow 
channel  configuration.  A  heat  transfer  model  coupling  conduction  through  the  thermoelement  legs  with 
convection  to  and  from  the  A1203  plates  is  formulated  to  investigate  the  performance  of  the  stack  as  func¬ 
tion  of  the  following  parameters:  hot  water  inlet  and  outlet  temperatures  (313-413  K  and  303-393  K, 
respectively),  stack  length  (300-1500  mm),  thermoelement  leg  length  (0.5-4  mm)  and  hot/cold  channel 
heights  (0.2-2  mm).  The  open-circuit  voltages  resulting  from  the  temperature  differences  are  within  3% 
mean  relative  error  of  those  resulting  from  temperature  differences  computed  via  CFD.  The  heat  transfer 
model  is  then  applied  to  optimize  a  1  kWe  stack  with  hot  water  inlet  and  outlet  temperatures  of  413  K 
and  393  K,  respectively,  for  either  a  maximum  heat-to-electricity  efficiency  of  4.2%  or  for  a  minimum 
volume  of  0.0021  m3. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

In  2010,  the  globally  installed  thermal  power  of  geothermal 
low-temperature  heat  (<150  °C  at  a  borehole  depth  of  <4  km) 
was  50.6  MWe  with  an  annual  growth  rate  of  12.3%  [1—3).  Geother¬ 
mal  low-temperature  heat  has  been  commonly  used  for  district 
heating  or  driving  heat  pumps.  As  of  recently,  it  has  also  been  con¬ 
sidered  for  electricity  generation  by  thermoelectric  converters 
(TEC)  [4-6], 
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A  TEC  module  comprises  p-type  and  n-type  semiconductor 
thermoelement  legs  fastened  between  two  ceramic  hot/cold  plates 
and  connected  in  parallel  thermally  and  in  series  electrically  [7-9], 
The  temperature  gradient  across  the  legs  induces  a  voltage  due  to 
the  Seebeck  effect  [10],  The  TEC  performance  is  characterized  by  its 
figure-of-merit  (ZT)  which  is  a  measure  of  the  heat-to-electricity 
efficiency:  the  higher  the  ZT,  the  higher  the  heat-to-electricity 
efficiency  is,  approaching  the  Carnot  efficiency  for  ZT  -»•  oo.  A 
state-of-the-art  BiTe  based  module  having  ZT  =  1  and  operating 
at  200  K  temperature  difference  reaches  the  theoretical  maximum 
efficiency  of  5%  [11], 

Thermoelectric  stacks  comprising  arrays  of  TEC  modules  adher¬ 
ing  to  parallel-,  counter-  and  cross-flow  water  heat  exchangers 
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Nomenclature 

A 

surface  (m2) 

Greek  letters 

B 

domain  boundary  (m2) 

oi,/S 

correction  factors 

Cp 

heat  capacity  (J  kg-1  K  ’) 

y 

iteration  step 

/ 

friction  factor 

error 

h 

heat  transfer  coefficient  (W  m  2  K_1) 

K 

thermal  conductivity  (W  m  1  K  ’) 

H 

height  (m) 

n 

heat-to-electricity  efficiency 

j 

electric  current  density  (A  m  2) 

p 

density  (kg  nr3) 

J 

electric  current  (A) 

p 

electrical  resistivity  (C2  m) 

l 

leg  length  (m) 

P 

electrical  contact  resistivity  (D  m2) 

L 

stack  length  (m) 

O 

power  density  (W  nr3) 

m 

water  mass  flow  (kg  s  ’) 

e 

pressure  (Pa) 

M 

number  of  thermoelement  legs 

n 

normal  vector 

Subscripts 

N 

number  of  modules 

c 

cold 

P 

electric  power  (W) 

conv 

convection 

Q 

heat  flux  (W  m  2) 

e 

electric 

a 

heat  transferred  (W) 

h 

hot 

R 

electrical  resistance  (fl) 

i 

finite  volume  index 

R2 

coefficient  of  determination 

leg 

thermoelement  leg 

S 

Seebeck  coefficient  (V K  ’) 

m 

middle 

t 

thickness  (m) 

OC 

open-circuit 

T 

temperature  (K) 

seg 

stack  segment 

u 

velocity  (ms1) 

w 

water 

U 

V 

voltage  (V) 
volume  (m3) 
volume  flow  (m3  s  ’) 
width  (m) 

Abbreviations 

V 

CFD 

computational  fluid  dynamics 

w 

ZJ 

TEC 

figure-of-merit 
thermoelectric  converter 

have  been  previously  analyzed  [5,12-16],  The  effect  of  the  dimen¬ 
sions  of  the  TEC  modules  and  flow  channels  on  the  conversion 
efficiency  has  been  examined  for  a  stack  with  counter-flow  config¬ 
uration  [5,16],  This  paper  expands  on  previous  studies  by  includ¬ 
ing  into  the  analysis  the  effect  of  the  geothermal  reservoir 
temperature,  i.e„  heat  transfer  fluid  inlet  temperature,  on  the  per¬ 
formance  of  a  1  kWe  stack  comprising  arrays  of  TEC  modules  with 
counter-flow  heat  exchangers.  A  heat  transfer  model  coupling 
conduction  through  the  thermoelement  legs  with  convection  to 
and  from  the  boundary  plates  is  formulated  after  evaluating  the 
predominant  modes  of  heat  transfer.  It  further  accounts  for 
temperature-dependent  physical  properties  and  heat  transfer 
coefficients  across  flow  regimes.  The  model  was  numerically 
solved  to  analyze  the  heat-to-electricity  efficiency  >;sta Ck  and  stack 
volume  \4tack  as  a  function  of  a  wide  hot  water  inlet  and  outlet 
temperature  range  and  geometrical  parameters.  Optimization  of 
the  stack  configuration  is  accomplished  for  either  (a)  maximum 
Vstack  or  (b)  minimum  r)stack  (maximum  power  density). 

2.  Scheme  of  thermoelectric  stack 

The  modular  stack  configuration  is  schematically  shown  in 
Fig.  1.  The  stack  comprises  lateral  arrays  of  TEC  modules.  Each 
module  is  represented  by  a  pair  of  A1203  plates  fastening  127  pairs 
of  p/n- type  (Bio.2Sbo.8)2Te3/Bi2(Teo.96Seo.o4)3  thermoelement  legs 
(M  =  254)  [17],  The  surface  of  the  A1203  plates  and  the  cross  section 
of  the  thermoelement  legs  are  Apiate  =  30  x  30  mm2  and 
Aieg=  1.05  x  1.05  mm2,  respectively;  the  lengths  of  the  stack  (L) 
and  the  thermoelement  legs  (1)  are  allowed  to  vary  as  model 
parameters.  A  temperature  gradient  across  the  thermoelement  legs 
is  imposed  via  a  parallel-plate  heat  exchanger  separating  the 
thermoelectric  arrays  and  operating  hot  and  cold  water  in  the 
counter-flow  configuration.  The  heat  exchanger  is  made  of  1  mm 
thick  Cu  plates  adhering  to  the  0.5  mm  thick  A1203  plates. 


3.  Heat  transfer  model 

The  heat  transfer  model  has  been  developed  for  a  stack  segment 
outlined  by  the  dashed  line  in  Fig.  1 .  A  more  detailed  schematic  of 
the  segment  is  shown  by  Fig.  2  representing  it  as  an  array  of  Nlrray 
TEC  modules  having  the  total  length  I,  the  width  W=  30  mm,  and 
hot  and  cold  water  half-channels  with  thickness  th  and  tc,  respec¬ 
tively,  resulting  in  the  total  height  H  =  l  +  th  +  tc  +  2-  tpiate  +  2  ■  tCu. 
The  following  assumptions  have  been  made:  (i)  heat  transfer  by 
radiation  and  free  convection  between  hot  and  cold  A1203  plates 
is  neglected  as  their  contribution  compared  to  the  total  heat  trans¬ 
fer  has  been  estimated  by  CFD  to  be  less  than  5%  (Appendix  A),  (ii) 
the  thermal  resistances  for  the  heat  conduction  through  the  Cu  and 
A1203  plates  is  neglected,  and  (iii)  the  hot  and  cold  sides  of  the  legs 
within  a  module  have  uniform  temperatures.  Thus,  the  heat  trans¬ 
fer  is  considered  to  occur  via:  (1)  convection  between  water  and 
the  thermoelement  legs,  and  (2)  ID  conduction  through  the  ther¬ 
moelement  legs.  Furthermore,  it  is  assumed  that  the  hot  and  cold 
sides  of  the  legs  within  a  module  have  uniform  temperatures. 
The  model  only  considers  steady-state. 

3.1.  Governing  equations 

Heat  transfer  is  formulated  for  the  hot  side  only  due  to  the  com¬ 
plete  analogy  with  the  cold  side.  The  general  continuous  energy 
conservation  equation  in  integral  form  for  the  hot  water  channel 
considering  heat  transport  by  the  water  flow  and  heat  convection 
through  walls  is  given  in  terms  of  the  hot  water  temperature 
TWih,  the  velocity  vector  of  the  water  flow  Uh,  the  normal  vector 
of  the  channel  domain  nh,  the  convective  heat  flux  gconvh  and  the 
channel  domain  boundary  B  by  [18]: 

Jb  Pvj  ■  cP,w  •  fw,h  •  (th,  ■  Uh)  •  dB  qconVih  nh  dB 


(1) 
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Fig.  1.  Schematic 
stack  segment  fo 


thermoelectric  stack  with  arrays  of  TEC  modules  separated  by  Cu  plate  heat  exchangers  in  counter-flow  configuration. ' 
l  the  heat  transfer  model  has  been  developed. 
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Fig.  2.  Schematic  of  the  modeled  stack  segment:  an  array  of  stack  length  I  with  Narray  TEC  modules  having  a  leg  length  1  between  hot  and  cold  water  half-channels  of 
thickness  th  and  tc\  respectively,  resulting  in  the  height  H.  The  channels  are  discretized  into  Narlay  cells  (dotted  boxes),  having  the  uniform  water  temperatures  Twh.i  and  TWpCii 
within  the  cell  i,  respectively.  The  dashed-dotted  lines  indicate  symmetry  boundaries. 


By  applying  the  finite  volume  method,  the  hot  water  channel  is 
discretized  into  Nar ray  rectangular  cells  (control  volumina,  see 
dotted  boxes  in  Fig.  2),  having  the  uniform  discrete  temperature 
Tw,h,i  within  a  single  cell.  For  the  discretization  the  following 
relations  are  used:  mseg,h  =  pw  ■  Uh  ■  ^channel,  ^channel  =  W  ■  th  and 
<Jconv.h  =  hh.i  ■  (Tw,h,i  -  Tieg  h  i).  Then  the  discretized  equation  in  the 
cell  i  in  terms  of  hot  water  mass  flow  mseglh,  hot  water  temperature 
Tw,h,i  and  hot  side  leg  temperature  T|eghi  yields  [19]: 

/Twhj_i  +  Twhi  Twhi  +  Twhi+1\ 
cP,w,,mseg,h  "  a"**) 

=  hh,i  -  Aplate  -  (Tw,h,i  -  Tleg,h, i)  (2) 

The  convective  (hot  water)-to-(thermoelement  leg)  heat  transfer 
coefficients  f!h,i  have  been  determined  from  Nusselt  correlations 
for  laminar,  transitional  and  turbulent  flows  between  parallel  flat 
plates  [20-23],  The  continuous  energy  conservation  equation  in  dif¬ 
ferential  form  for  the  thermoelement  legs  considering  Fourier’s  heat 
conduction,  Joule  heating,  and  Thomson  effect  is  given  in  terms  of 
the  temperature  in  the  leg  Tleg  and  the  current  density  j  by: 

v(/cleg  ■  VTleg)  +  pleg  f  -  Tteg  Ses  j-  VTleg  =  0  (3) 

ai  leg 

where  the  current  density  j  satisfies  current  conservation  Vj  =  0 
and  is  given  in  terms  of  the  voltage  Uleg  and  the  leg  temperature  Tieg 
by: 


[VLfieg  +  Sleg  ■  Vr,eg] 


(4) 


With  the  assumption  of  1 D  heat  conduction  in  the  thermoelement 
legs,  linear  dependency  of  the  thermoelectric  material  properties 
on  temperature  and  the  overall  current  J,  Eq.  (3)  can  be  algebraically 
solved.  Thus,  the  ID  heat  conduction  at  the  hot  side  of  a  thermoel¬ 
ement  leg  in  the  cell  i  in  terms  of  the  hot  side  leg  temperature  Tiegih,i, 
the  overall  current  J  and  the  temperature  difference  across  the  legs 
(ATieg,i  =  Tieg,h,i  -  Tieg,c,i)  reduces  to: 


Qleg,W  =  S,eg,m,i  ■  ^eg.h.i  J ~  «h,i '  l  ^  f/ 2  +  A,,i '  Aeg^‘  ATlegii 


where 

«h,i=  (1- 


1  ^Pleg.m.i  ATiegA 
P,eg.m,i  '  dhee,m:i  '  4  J 


are  the  correction  factors  accounting  for  the  linear  dependency  of 
the  material  properties  on  temperature  [24],  As  aforementioned, 
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the  heat  transfer  considerations  for  the  cold  side  are  analogous  to 
the  ones  presented  above.  Further  the  electric  current  is  given  by: 


J  =  ff array  /  (Rarray  +  ^contact  +  Rload)  (8) 

The  voltage  is: 

Uarray  =  M.gVSleg,m,i-Arieg,i  (9) 

the  material  resistance: 

Rarray  =  M.^£VAeg.m,i  (10) 

and  the  contact  resistance: 

Rcontact  =  2  Narray  M-/-  (11) 

/ileg 

with  the  matched  load: 

Rload  =  Rcontact  +  Rarray  (12) 

The  power  output  of  the  TEC  module  array  is  calculated  by: 

Rarray  =/-Rload  (13) 


Physical  properties  of  water  and  thermoelectric  properties 
of  the  thermoelements  are  estimated  by  polynomial  correlations 
[25], 

3.2.  Boundary  conditions 

The  power  output  of  the  stack  was  set  to  Pstack  =  1  kWe.  The  hot 
water  inlet  (TWihiin)  and  outlet  (TWih  out)  temperatures  were  varied  in 
the  range  TWih,in  =  313-413  K  and  fw,hiOUt  =  303-393  K.  respectively, 
yielding  hot  water  temperature  differences  (ATw,h)  in  the  range 
ATWh  =  10-110  K.  The  cold  water  inlet  (TWiC,in)  and  outlet  (Tw,c, out) 
temperatures  were  set  to  Twcin  =  293  K  and  TWiCi0ut  =  298  K,  respec¬ 
tively,  yielding  the  cold  water  temperature  difference  ATWC  =  5  K. 
The  mass  flow  rates  of  hot  (mstack,h)  and  cold  water  (mstack,c)  were 
adjusted  to  satisfy  the  prescribed  hot/cold  water  outlet  tempera¬ 
tures.  The  axis  of  symmetry  was  set  at  the  middle  of  the  water 
channels  (dashed-dotted  line  in  Fig.  2).  The  pressure  at  the  inlets 
to  the  channels  is  5  bars. 

3.3.  Model  parameters 

The  geometric  parameters  were  ranging  as  follows: 


and  the  useful  power  output  of  the  stack  segment  is  given  by: 

Pseg  =  Rarray  -  Ppump  (14) 

where  the  pump  work  RpUrnp  is  estimated  in  terms  of  the  pressure 
drop  in  the  water  channels  and  the  volume  flow.  The  pressure 
drop  in  the  hot  channel  -  analogously  for  the  cold  channel  -  is 
estimated  in  terms  of  the  friction  factor  /h  and  the  flow  velocity 
Uh  by: 

A©h=/h^  (~-pw-Uh)  05) 

The  friction  factor  fh  has  been  determined  for  laminar,  transitional, 
and  turbulent  flows  between  parallel  flat  plates  [20,22],  Then  the 
pumping  work  in  the  hot  channel  (analogously  for  the  cold  channel) 
is  calculated  in  terms  of  the  pressure  drop  A0h  and  the  volume 
flow  Vh  by: 


(1 )  stack  segment  length  L  =  300-1500  mm,  i.e.,  Narray  =  10-50 

(2)  leg  length  /  =  0.5-4  mm  and 

(3)  hot/cold  half-channel  thicknesses  {th,tc}  =  0.1-1  mm. 

The  model  parameters  are  shown  in  Table  1 . 

3.4.  Numerical  solution 

The  coupled  governing  Eqs.  (1)— (19)  were  implemented  in  a 
FORTRAN  code  and  solved  iteratively  until  the  convergence  crite¬ 
rion  |  (Tj,  h  j  -  j) /Tj, h  i|  g:  10  6  for  all  elements  (i)  and  the  overall 
energy  balance  |(AQh  -  AQC  -  Pseg)/AQh|  ^  1(T4  were  satisfied, 
where  AQh  =  mseg,h  ■  cp,w  ■  ATw,h  and  ACL  =  mseg,c  ■  cp,w  ■  ATWiC. 

3.5.  Comparison  with  CFD-based  model 


Rpump.h  —  '  A©h 


(16) 


Evidently,  the  total  pumping  work  is  then  Rpump  =  Rpunlp,h  +  Rpump.c- 
Finally,  the  heat-to-electricity  efficiency  of  the  stack  segment  tjseg 
is  defined  as  the  ratio  of  the  electrical  power  output  over  the  heat 
delivered  to  the  stack  segment: 


Pse  g 

tpft  \r»6; 


(17) 


The  volume  of  the  stack  segment  is  given  in  terms  of  the  stack 
length  L,  the  stack  segment  height  H,  and  the  stack  segment  width 
Why: 


Vseg  =  L-H-W 


(18) 


The  power  density  of  the  stack  segment  is  expressed  as: 


(19) 


Finally  for  the  total  stack,  the  magnitudes  t/sta ck.  Vstack,  mh  stack  Nstack, 
and  o stack  are  calculated  in  terms  of  the  ratio  of  the  power  outputs 
of  the  stack  Pstack  and  stack  segment  Pseg  by: 


Nstack  =  »7seg.  ^ stack  =  VSe| 


=  mhseg  ■  paCk  ,  Nstack=Narr. 


(20) 


The  steady-state  momentum  and  heat  balance  equations  for  the 
stack  segment  have  been  solved  using  computational  fluid 


rwt,,in  313-413 

rw'h'out  303-393 

Ae g  1-05  x  1.05 

Apiate  30  X  30 

l  0.5-4 

L  0.3-1. 5 

Narray  10-50 

M  254 

tplate  0.5 

tcu  1 

tc  0.1-1 

th  0.1-1 

Kp-type  1.13-1.36 

Sp-type  204-235 

Pp-,ype  10.4-16.3 

K„_type  1.26-1.38 

S„-,ype  -215  to  -207 

Pn-type  11-15.5 

p  5e-9 


Unit 


kW 

K 

K 


W  m~’ 
|iVK  1 

pn  mi 


pVK-* 1 

pnm 

nm2 


Source 

Set 

Model  parameter 
Model  parameter 
Assumed 
Assumed 
Model  parameter 
Model  parameter 
Model  parameter 
Set 

Assumed 
Model  parameter 
Model  parameter 
Set 

Assumed  [17] 
Assumed  [17] 
Assumed  [17] 
Assumed  [17] 
Assumed  [17] 
Assumed  [17] 
Assumed  [27] 
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R2  =  0.997.  The  mean  relative  error  is  £  =  3%. 


rw,h,in[K] 


function  of  the  hot  water  inlet  and  outlet  temperature  for  the  case  of  maximized 
heat-to-electricity  efficiency. 


dynamics  (CFD)  software  ANSYS  CFX  13.0  [26]  in  order  to  obtain 
the  rigorous  temperature  profiles  of  cold  and  hot  sides  of  the  legs 
along  a  segment.  Open  circuit  voltages  were  then  calculated  using 
the  mean  temperature  differences  across  the  TEC  modules  and 
compared  in  Fig.  3  with  those  calculated  based  on  the  heat  transfer 
model  described  in  this  work.  The  high  coefficient  of  determination 


Fig.  5.  (a)  The  stack  volume  and  the  power  density  as  function  of  the  hot  water  inlet 
and  outlet  temperature;  and  (b)  the  stack  volume  as  function  of  the  number  of  TEC 
modules  in  the  stack  for  the  case  of  minimized  stack  volume. 


in  this  parity  plot  ( R 2  =  0.997)  implies  good  agreement  between  the 
two  methods,  further  supported  with  mean  relative  error  of  e  =  3%. 

4.  Simulations  of  a  1  lcWei  stack 

A  stack  delivering  the  useful  power  output  Pstack  =  1  kWei  has 
been  optimized  by  varying  geometric  parameters  at  a  given  set 
of  hot  water  inlet  (TWih,in)  and  outlet  (TWih,out)  temperatures  having 
either  (a)  maximum  heat-to-electricity  efficiency  or  (b)  minimum 
size  (maximum  power  density)  as  goal  function. 

4.1.  Maximum  heat-to-electricity  efficiency 

The  heat-to-electricity  efficiency  (f? stack)  as  function  of  TWih,,n 
and  rWih,out  is  shown  in  Fig.  4a.  It  almost  linearly  increases  with 
Tw,h,in  and  Tw,h,out.  implying  a  stronger  dependency  on  the  latter. 
The  calculated  values  range  from  f/sta Ck  =  0.5%,  obtained  for 
Tw,h.in  =  313  K  and  Tw,h>out  =  303  K;  to  rjs tack  =  4.2%  for  rw.h,in  =  413  K 
and  Tw  h  out  =  393  K. 

Fig.  4b  presents  the  effect  of  TWih  in  and  TWihiOUt  on  the  hot  water 
mass  flow  (mstack,h).  mstack,h  monotonically  decreases  with  TWih,in  but 
increases  with  TWih,0ut-  The  values  range  from  mstack,h  =  4.85  kg/s, 
calculated  for  Tw,h,in  =  313  K  and  Tw,h,out  =  303  K,  to  mstack>h  =  0.28 
kg/s,  calculated  with  TWihiin  =  413  K  and  TWih,out  =  303  K.  One  should 
notice  that  increase  in  TWihJn  leads  to  an  increase  in  f/stack  but  a 
decrease  in  mstackih  while  both  jjstack  and  mstack,h  increase  with  an 
increase  in  Tw  h,0ut. 
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Table  2 


(a)  Maximized  r/sack  (b)  Minimized  Vstack  Comparing  (a)  to  (b) 

Jw.h.o„,  T'w.h.in  'istack  ^stack  <Wk  to  ^stack  <Tstack  lis,ack.(a)/'7s,ack.(b)  Vstack.(a)/ l/stack.(b)  (<X5tack.(b)/ to.(a)) 

(K)(K)  (%)  (m3)  (kWm-3)  (%)  (m3)  (kWm-3)  (-) (-) 

313  303  0.5  0.9075  1.1  0.16  0.2871  3.5  3.13  3.16 

413  303  1.9  0.0611  16.4  0.78  0.0146  68.5  2.44  4.18 

413  393  4.2  0.0135  74.1  2.1  0.0021  476.0  2.00  6.41 


4.2.  Minimum  stack  volume 

The  stack  volume  (V4tack)  and  the  power  density  (<7stack)  ate 
shown  in  Fig.  5a  as  function  of  TWih.,n  and  Tw  h  out-  Vstack  monotoni- 
cally  decreases  and  attack  monotonically  increases  with  Tw,h,in  and 
fw.h.out-  The  stack  volumes  range  from  VstaCk  =  0.29  m3  (<7stack  =  3.5 
kW  m  3),  calculated  for  Tw>h,in  =  313K  and  Tw>h,out  =  303  K,  to 
Vstack  =  0.0021  m3  (CTstack  =  476kWm  3)  obtained  with  TWih,in  = 
413  K  and  Tw,h,out  =  393  K. 

Fig.  5b  is  showing  \4tack  as  function  of  the  number  of  modules 
(Nstack)  for  all  simulated  hot  inlet  and  outlet  temperatures.  The 
data  points  indicate  a  linear  relationship  (correlation  factor  of 
0.9947).  The  number  of  modules  ranges  from  JVstac k  =  550  to  JVstac k  = 
79,740  as  calculated  with  l/stack  =  0.0021  m3  and  VstaCk  =  0.29  m3, 
respectively. 


5.  Summary  and  conclusions 

A  heat  transfer  model  coupling  ID  conduction  through  the  ther¬ 
moelement  legs  with  convection  to  and  from  the  legs  has  been  imple¬ 
mented  in  order  to  investigate  the  performance  of  a  thermoelectric 
stack  as  function  of  operating  parameters  and  stack  geometry.  The 
simulated  open-circuit  voltages  have  been  found  to  be  within  3% 
mean  relative  error  as  compared  to  those  calculated  via  CFD. 

A  stack  delivering  the  useful  power  output  Pstack  =  1  kWe  has 
been  optimized  by  varying  geometric  parameters  at  a  given  set 
of  hot  water  inlet  (TWih,i„)  and  outlet  (TWih,out)  temperatures  having 
as  the  goal  function  either  (a)  maximum  efficiency  or  (b)  minimum 
size  (maximum  power  density).  The  key  results  are  summarized  in 
Table  2,  implying  strong  dependency  of  both  t]stack  and  Vstack 
(ffstack)  on  Twh.in  and  Tw,h,out. 

The  results  have  opened  the  possibility  to  perform  economic 
analyses  considering  the  costs  of  the  TEC  stack  and  the  geothermal 


system  in  order  to  define  the  configuration  and  operating  parame¬ 
ter  range  that  are  optimal  from  a  commercial  standpoint. 

Appendix  A 

The  magnitude  of  convective  and  radiative  heat  exchange  be¬ 
tween  the  hot  and  cold  plates  relative  to  the  heat  conducted 
through  the  thermoelement  legs  has  been  evaluated  by  CFD  based 
on  the  following  worst-case-scenario,  chosen  to  overstate  the  con¬ 
tributions  of  convective  and  radiative  transfers: 

(1 )  Hot  side  temperature  of  413  K  (the  maximum  considered  by 
simulations). 

(2)  Cold  side  temperature  of  293  K  (the  minimum  considered  by 
simulations). 

(3)  Leg  length  of  4  mm  (the  longest  considered  by  simulations). 

The  CFD  model  considers  a  4-leg  module.  2  p/n- type  thermoele¬ 
ment  leg  pairs  and  the  ambient  air  are  simulated,  as  shown  in 
Fig.  Al.  The  thermoelement  legs  have  dimensions  of  4  x  1.05  x 
1.05  mm3  and  a  distance  of  0.92  mm  between  neighboring  legs. 
Both,  the  dimensions  of  the  legs  and  the  distance  between  neigh¬ 
boring  legs  correspond  to  those  for  the  investigated  254-leg 
modules.  The  thermal  resistance  of  the  A1203  plates  is  neglected, 
i.e.  they  have  no  thickness,  and  constant  temperature  conditions 
are  imposed  on  the  hot  (Th  =  413  K)  and  cold  side  (Tc  =  293  K)  of 
the  thermoelement  legs  (indicated  by  solid  lines).  Symmetry 
boundaries  are  set  at  the  lateral  sides  (indicated  by  dashed-dotted 
lines).  Radiative  heat  exchange  is  simulated  by  the  Monte  Carlo 
method  using  109  rays.  The  emissivities  of  the  A1203  plates  and 
thermoelement  legs  are  0.3  and  0.7,  respectively.  The  CFD  mesh 
contains  291,254  elements  where  thep/n-type  thermoelement  legs 
and  the  air  domain  have  hexahedral  and  tetrahedral  cells,  respec¬ 
tively.  Convergence  was  achieved  when  the  residuals  were  <10  6. 
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